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Experimental signals of non-linear magneto-optical resonances at Di excitation of natural ru- 
bidium in a vapor cell have been obtained and described with experimental accuracy by a detailed 
theoretical model based on the optical Bloch equations. The Di transition of rubidium is a challeng- 
ing system to analyze theoretically because it contains transitions that are only partially resolved 
under Doppler broadening. The theoretical model took into account all nearby transitions, the 
coherence properties of the exciting laser radiation, and the mixing of magnetic sublevels in an 
external magnetic field and also included averaging over the Doppler profile. Great care was taken 
to obtain accurate experimental signals and avoid systematic errors. The experimental signals were 
reproduced very well at each hyperfine transition and over a wide range of laser power densities, 
beam diameters, and laser detunings from the exact transition frequency. The bright resonance 



expected at the Fg — 1 
observed at the F„ — 2 



> Fe — 2 transition of ^'^Rb has been observed. A bright resonance was 
Fe — i transition of *^Rb, but displaced from the exact position of the 
transition due to the influence of the nearby Fg = 2 Fe = 2 transition, which is a dark resonance 
whose contrast is almost two orders of magnitude larger than the contrast of the bright resonance 
at the = 2 — > _Fe = 3 transition. Even in this very delicate situation, the theoretical model 
described in detail the experimental signals at different laser detunings. 

PACS numbers: 32.60.-(-i,32.80.Xx,32.10.Fn 
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I. INTRODUCTION 



Dark [l|, Q and bright [1] non-linear magneto-optical 
resonances associated with atomic state coherences 
and strongly related to the ground state Hanle effect [5] 
(see |6|] for a review) are a particular manifestation of 
a broad class of non-linear magneto-optical phenom- 
ena (see, for example, 0, 's'l for a review). Dark reso- 
nances in a scan of laser induced fluorescence (LIF) ver- 
sus magnetic field arise when at zero magnetic field co- 
herences are induced among ground state magnetic sub- 
levels and a quantum superposition state is created that 
is not coupled to the exciting laser field. Bright reso- 
nances arise when atoms at zero magnetic field are pre- 
pared in a superposition of ground state sublevels that 
is strongly coupled to the excited state. These reso- 
nances have a variety of applications, such as in op- 
tical magnetometers Q and optical switches and 
they have been studied extensively in many systems 
both theoretically and experimentally (see, for exam- 
ple, |,|l,|illil,[il,|ii,[il,|li[i3,lillli,H|2li). How- 
ever, accurately modelling the profiles of bright and dark 
resonances in systems that are unresolved under Doppler 
broadening remains a challenging enterprise. There are 
two ways to address the issue of Doppler broadening in 
a detailed study of these resonances in alkali vapors: (1) 
compensating for the lack of spectral resolution by us- 
ing a detailed theoretical model or (2) having recourse 
to an extremely thin cell whose width is comparable to 



the wavelength of the laser radiation, and which has sub- 
Doppler resolution [2l[ . We chose the first approach and 
aimed to demonstrate that it is possible to achieve with 
experimental accuracy a quantitative description of the 
shape of all bright and dark resonances of the Di tran- 
sition of ^^Rb and ^^Rb by means of a straightforward 
theoretical model based on the optical Bloch equations 
with minimal recourse to adjustable parameters. In par- 
ticular, we wanted to study the dependence of the res- 
onance shapes, in particular the resonance contrast and 
width, on the laser power density, laser frequency, and 
laser beam diameter. The laser beam diameter is related 
to the transit relaxation time, which is critical for de- 
termining the resonance profile. Previous works did not 
study how the resonance shape varies with beam diame- 
ter. 

The ground state Hanle effect was first observed by 
Lehmann and Cohen- Tannoudji in 1964 Q. In 1978, 
Pique presented measurements of dark resonances in the 
Fg = 2 — > i^e = 1 transition of the Di line in a beam of 
sodium atoms with negligible Doppler broadening, which 
agreed quite well with calculations based on the opti- 
cal Bloch equations 11 1. Renzoni ct al. [l^] studied ex- 
perimentally and theoretically all transitions of the Di 
sodium line in a beam of sodium atoms, but saw only 
dark resonances. Bright resonances were first observed 
in 2000 in the Di and D2 transitions of rubidium by 
Dancheva et al. i3| in a vapor of rubidium atoms, but 
the study did not include a model-based analysis of the 
results. Papoyan et al. jlH studied bright and dark res- 
onances in the D2 transition of cesium, but their theo- 
retical model based on the optical Bloch equations pre- 
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transition, which appeared to be dark in the experiment. 
Alzetta et al. [l^ excited the Di transition of sodium 
in a vapor cell with a broadband multi-mode laser to 
demonstrate dark resonances with 100% contrast and 
modelled these resonances with rate equations. However, 
the model was not used to predict the shape of the res- 
onances. Gateva et al. measured the widths of the 
dark resonances at the Fg = 2 ^ Fe = 1 transition of the 
Di line of ^^Rb and compared the widths to predictions 
based on an analytical expression obtained for a Doppler 
broadened three-level system in the linear approxima- 
tion. Most recently (2ll | , experiments were performed on 
the D2 transition of cesium in an extremely thin cell, 
but some resonances that should have been bright ap- 
peared dark for as yet unknown reasons, and the reso- 
nances are rather broad. The most advanced previous 
attempt to apply a detailed model magneto-optical reso- 
nances in alkali vapors was made by Andreeva et al. [l^ , 
who worked on the cesium D2 transition, which is un- 
resolved under Doppler broadening. Their model took 
into account all nearby hyperfine transitions and aver- 
aged over the Doppler profile and over the distribution 
of transit times through the laser beam. Nevertheless, 
the theoretical and experimental contrast of these reso- 
nances as a function of laser power density presented in 
this study differed by about a factor of two. Moreover, 
the paper presented a comparison of theoretical and mea- 
sured resonance contrast as a function of laser power den- 
sity only for one bright resonance and did not study the 
dependence on beam diameter. One weak point in the 
model was that it did not take into account the coherence 
properties of the laser radiation in sufficient detail. 

In this paper we endeavored to show that the optical 
Bloch equations can model correctly bright and dark res- 
onances in vapor cells when the model takes into account 
nearby hyperfine levels and magnetic sublevel mixing, av- 
erages over the Doppler profile, and carefully includes the 
coherence properties of the laser radiation. We first had 
applied such a detailed model to the Di line of cesium in 
an atomic vapor [2^ . Cesium Di excitation has the sim- 
plifying advantage that the energy splitting between the 
excited state hyperfine levels exceeds the Doppler broad- 
ening. After having gained confidence in the model's suit- 
ability for describing a simple situation, we now report 
on studies of a more complex system, which tested more 
rigorously the applicability of the model under a wide 
variety of experimental conditions. We chose the Di line 
of natural rubidium, whose excited state hyperfine levels 
can be resolved only partially under Doppler broadening. 
Figuredlshows the level scheme of the Di line in ^^Rb and 
®^Rb. The energy difference between the excited state 
hyperfine levels is 361.6 MHz for ^^Rb and 816.7 Mhz for 
s^Rb. The full width at half maximum (FWHM) of the 
Doppler broadening for rubidium at room temperature is 
about 500 MHz. The arrows represent the possible tran- 
sitions between ground and excited states, and the frac- 
tions next to the arrows represent the relative transition 
strengths. Figure [5] shows a LIF excitation spectrum of 



the Di line of natural rubidium obtained with the same 
experimental equipment as the other results described in 
this paper. The exact positions of the resonances can be 
seen from the simultaneously recorded saturated absorp- 
tion spectrum. One can see that the two excited state hy- 
perfine levels of *^Rb are almost resolved under Doppler 
broadening, while the hyperfine levels of ^^Rb are rather 
poorly resolved, which is consistent with the energy level 
separations given in Figure [T] and the Doppler width of 
rubidium at room temperature. This system is interest- 
ing, because it is expected to contain six dark and two 
bright resonances, of which the bright resonance at the 
Fg = 1 ^ Fe = 2 oi ^^Rb has never before been observed 
as far as we know. 
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FIG. 1: Hyperfine level structure and transitions of the Di 
line of rubidium. The fractions on the arrows indicate the 
relative transition strengths [23I . [24l ] . 
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FIG. 2: Spectrum of Rb Di line. Curve 1 — saturated ab- 
sorption spectrum; curve 2 — LIF excitation spectrum at 10 
mW/cm^ excitation. 



II. EXPERIMENTAL DESCRIPTION 

The experimental geometry is shown in Figure[31 which 
depicts the relative orientation of the laser beam (exc), 
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the laser radiation's linear polarization vector (Eea;c)7 the 
magnetic field (B), and the observation direction (obs). 
The experimental setup and auxiliary equipment are sim- 
ilar to what was used in [22*1. In the present experiment, 
natural rubidium was confined at room temperature in 
cylindrical glass cells. The laser was a DL 100 exter- 
nal cavity single mode diode laser with a wavelength of 
794.3 nm and a typical line width of a few Megahertz, 
which was produced by Toptica, A. G., of Graefelfing, 
Germany. Two different cells were used in the experi- 
ments. One was produced by Toptica, A. G., and made 
of Pyrex with optical quality windows with a diameter of 
25 mm. The other cell was made in our laboratory and 
had a diameter of 35 mm. The cell was placed at the 
center of a three-axis Helmholtz coil system. Two pairs 
of coils compensated the laboratory magnetic field, while 
the third pair of coils was used to scan the magnetic field 
in the observation direction. The current in this third 
pair of coils was scanned with a Kepco BOP-50-8M bipo- 
lar power supply that was controlled by an analog signal 
from a computer. For some of the experiments, part of 
the laser radiation was diverted to the second cell, which 
was used to obtain a saturated absorption spectrum, from 
which the exact position of the transitions were deter- 
mined. In addition to the National Instruments 6024E 
data acquisition card, an Agilent DSO5014A oscilloscope 
was used in some of the experiments for data acquisition. 
Both data acquisition systems produced similar results. 

The laser induced fluorescence was detected by means 
of a Thorlabs FDS-100 photodiode as a function of mag- 
netic field for all hyperfine transitions of the two iso- 
topes of rubidium (see Fig. [T|) . No polarizers were used 
in the LIF observation. Signals were obtained for differ- 
ent laser power densities (between 5 mW/cm^ and 1000 
mW/cm'^) and laser beam diameters (between 0.09 mm 
and 2.3 mm). The beam diameter was assumed to be 
the full width at half maximum (FWHM) of the beam 
profile, which was determined by means of a Thorlabs 
BP 104- VIS beam profiler. The beam profile was roughly 
a symmetrical Gaussian. 

Signals were obtained by employing a double scanning 
technique [17i | in which the laser frequency was scanned 
slowly across a transition or pair of adjacent transitions 
while the magnetic field was scanned more rapidly follow- 
ing a saw-tooth shape from negative to positive values. 
The laser frequency changed by about 10-20 MHz per 
second. The typical laser frequency scan lasted about 
2-3 minutes, whereas the typical period for the magnetic 
field scan was 1 second. In this manner, a series of reso- 
nance signals at laser frequencies differing by 10-20 MHz 
could be obtained. For some observations the laser's fre- 
quency was maintained at the maximum of the fluores- 
cence signal at a given transition while the magnetic field 
was scanned several times, and the scans were averaged. 
The laser was not actively stabilized, but its frequency 
was monitored with a High Finesse WS-7 wavemeter and 
was found to vary by a quantity on the order of 10 MHz 
during the time of a typical measurement. 



In addition to the signal, each measurement included 
a certain amount of background. Background from scat- 
tered laser light could be determined easily by tuning the 
laser off resonance. However, the background associated 
with scattered LIF could not be readily identified. It 
was assumed that this background accounted for a fixed 
percentage of the signal for a given vapor cell. The ex- 
perimental signals from the home-made cell matched the 
calculated curves when the signals from the homemade 
cell were assumed to contain about 20% background from 
scattered LIF. Because of the improved transparency of 
the high optical quality glass in the cell from Toptica, 
A. G., the component of the background from scattered 
LIF was so small that the theoretical and experimental 
curves agreed without assuming any contamination from 
scattered LIF. With this cell, it was no longer necessary 
to consider the background from scattered LIF to be an 
adjustable parameter in the theoretical treatment of the 
signals. Almost all the signals were remeasured with the 
cell from Toptica, A.G., in order to verify our results. 
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FIG. 3: Experimental geometry. The relative orientation of 
the laser beam (esc), laser light polarization (Eexc), magnetic 
field (B), and observation direction (obs) are shown. 



III. THEORETICAL MODEL 

The theoretical model is basically the same as the one 
described in [13]. The model assumes that the atoms 
move classically and are excited at the internal transi- 
tions. Thus the internal atomic dynamics can be de- 
scribed by semi-classical atomic density matrix p, which 
also depends parametrically on the classical coordinates 
of the atomic center of mass. The time evolution of the 
density matrix p is described by the optical Bloch equa- 
tions (QBE) [ll 



H,p 



ihRp, 



(1) 



where the relaxation operator R includes the sponta- 
neous emission rate, which equals the natural transi- 
tion linewidth F and the transit relaxation rate 7 (as 
the atom-laser interaction takes place in a finite spatial 
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region). The Hamiltonian H can be written as follows 
H = Ho + Hb + V, where Hq is the unperturbed atomic 
Hamiltonian, which depends on the internal atomic co- 
ordinates, Hb is the Hamiltonian of the atomic interac- 
tion with the magnetic field, and V = — d • E{t) is the 
dipole interaction operator, which includes the electric 
dipolc operator d and the electric field of the excitation 
light E(t). The electric field is described classically with 
some polarization in the direction of the unit vector e, an 
electric field amplitude e^, a central frequency uJ, and a 
fluctuating phase ^(t), which causes a finite bandwidth 
Aoj (FWHM) for a Lorcntzian shaped spectral profile: 

E(t) = e{t)e + e*{t)e*, 



(2) 

e{t) = |e-|e-^*(*)-*(--k^-)*. (3) 

The atoms move with definite velocity v, which causes 
the shift kjjv in the laser frequency due to the Doppler 
effect, with k— being the wave vector of the excitation 
light. 

Then we apply the rotating wave approximation [26] to 
the OBEs. ft results in stochastic differential equations, 
which can be further simplified by using the decorrela- 
tion approach [27]. The stochasticity is caused by the 
random fluctuations of the laser radiation and leads to 
a finite spectral width of the radiation. In the experi- 
ment we observed the light intensity, which is a quan- 
tity that is averaged over time intervals that are large 
in comparison with the characteristic time of the phase- 
fluctuations. Therefore we perform a statistical averag- 
ing of the stochastic differential equations. This is done 
by taking a formal statistical average over the fluctuat- 
ing phases and then using the decorrelation approxima- 
tion [2^, which is valid for either the random phase jump 
or the continuous random phase diffusion model. As a 
result, we eliminate the density matrix elements that cor- 
respond to optical coherences and arrive at equations for 
the Zeeman coherences, 



dp. 



dt 
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where Pg^gj and pe^cj are the density matrix elements 
for the ground and excited states, respectively. The first 
term in ([4]) describes the re-population of the ground 
state and the creation of Zeeman coherences due to in- 
duced transitions, Fp^g^e^ and ^p^agj represent the laser 
field coupling of the ground and excited states, and d^'"^ 
is the dipole transition matrix element. The second term 
stands for the changes of ground state Zeeman sublevel 
population and the creation of ground state Zeeman co- 
herences due to light absorption. The third term de- 
scribes the destruction of the ground state Zeeman coher- 
ences by the external magnetic field; ojg^g^ is the splitting 
of the ground state Zeeman sublevels. The fourth term 
describes the re-population and transfer of excited state 
coherences to the ground state due to spontaneous tran- 
sitions. We assume our transition to be closed (within 
the hyperfine structure), and so ^ ^g'is] — T- The fifth 

ei,ej 

and sixth terms show the relaxation and re-population 
of the ground state due to non-optical reasons - in our 
case it is assumed to be solely transit relaxation. It is 
assumed that the atomic equilibrium density outside the 
interaction region is normalized to 1, therefore A = 7. 

In equation ^ the first term stands for the light ab- 
sorbing transitions; the second term denotes induced 
transitions to the ground state; the third describes the 
destruction of ground state Zeeman coherences in the 
external magnetic field; and the fourth term denotes the 
rate of spontaneous decay of the excited state, where 
LUeiej is the splitting of the excited state Zeeman sub- 
levels. 

The interaction strength Fp^g^e^ is given by 



I 



(6) 



where 



is proportional to the laser power density. 



The magnetic field not only leads to Zeeman sublevel 
splitting Wy, but changes the transition dipole elements 
due to magnetic sublevel mixing. In the case of two hy- 
perfine structure levels this mixing can be obtained from 
the Breit-Rabi formula 29, 30]. 

We are conducting experiments in stationary excita- 
tion conditions so that dpg^gjdt = dpeiejdt — 0. Thus, 
the differential equations -I- (O are reduced to a 
system of linear equations, which, when solved, yields 
the density matrices for the atomic ground and excited 
states. The observed fluorescence intensity can be calcu- 
lated as follows: 



^/(e)=/o 



Aob}*j(ob) 



(7) 



where Iq is a constant of proportionality. 

To include all the different velocity groups of atoms, we 
average over the Doppler profile simply by calculating the 
intensity for each velocity group and multiplying this in- 
tensity by the respective statistical weight corresponding 
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to the number of atoms in the respective velocity group. 
Also, to obtain the total signal for a specific observation 
direction, we average the fluorescence signal over the two 
orthogonal polarization components. Finally, since we do 
not resolve the fluorescence spectrally, we sum the signal 
over all four possible radiating hyperfine transitions. 

There are two parameters in the model which need 
to be determined by comparing the results from calcula- 
tions with experiment: (1) the conversion between laser 
power density / and Rabi frequency f2 and (2) the conver- 
sion between beam diameter and transit relaxation time 
7^^. Both of them can be theoretically estimated [sills^ . 
though the estimate of the conversion factor in the case of 
the Rabi frequency is rather rough and needed to be ad- 
justed by a factor of less than two to yield the best match 
between theory and experiment. The transit relaxation 
rate can be estimated from the mean thermal velocity of 
the atoms Vavg and the laser beam diameter d (FWHM) , 
and is close to Vavg/d. Since we assumed that fl cx \/T 
and 7 cx l/d, the same values for these two conversion 
factors were used in all our measurements, which, for two 
rubidium isotopes, all together spanned 8 hyperfine tran- 
sition and a large range of laser power densities and beam 
diameters. We also used 293 K for the temperature, 10 
MHz for the laser linewidth. Atomic constants such as 
nuclear spin Lande factor and the hyperfine structure 
constant A were taken from [H, [23| ■ Hyperfine structure 
constants B were assumed to be negligibly small. 



IV. RESULTS AND DISCUSSION 

Figures [3] and [5] show magneto-optical resonances at 
laser excitation frequencies tuned to the corresponding 
hyperfine transitions of the Di line of rubidium. Figured] 
shows the transitions of *^Rb at 100 mW/cm^, whereas 
Fig. [S] shows the transitions of ^^Rb at 10 mW/cm^. 
Experimental measurements are represented by markers, 
whereas the results of theoretical calculations are shown 
by solid lines. In each figure, panels (a)-(c) have the 
same vertical scale, but the vertical scale of panel (d) 
is much smaller. In most cases, the theoretical model 
reproduces the measured signals to within experimen- 
tal accuracy. Only at high laser power densities (> 100 
mW/cm^), some discrepancies arise as regards the reso- 
nance widths, although the model predicts the resonance 
contrast quite accurately even at high laser power den- 
sities. Indeed, in Figure UJ^a), which describes the ^^Rb 
Fg = 2 ^ Fe = 1 transition, the calculations predict a 
somewhat larger resonance width than is observed. This 
transition is the strongest transition at a rather high laser 
power density, and so absorption saturation effects could 
be strong enough at the wings of the spatial beam pro- 
file. At strong saturation of the absorption, it is no longer 
adequate to describe the transit relaxation rate by a sin- 
gle rate constant [12, H^, as is done in our theoretical 
model, since the resonance width is determined mainly 
by the ground state relaxation rate, which is dominated 



by the transit relaxation rate in this experiment. 

Figures |31Ia)-(c) and [5Ja)-(c) show dark resonances, 
which are expected, since in each case {Fg > Fe). The 
contrasts for the two cases where Fg = Fe (Fig. [D^b) and 
(c)) are equal, whereas when Fg > F^ (Fig. Sta)) the 
contrast is larger by more than a factor of two. Sim- 
ilar behavior is seen in Figure O although in this case 
the two transitions with Fg = F^ do not have similar 
contrasts, because adjacent resonances are not resolved. 
The i^g = 3 ^ Fe = 3 transition (Fig. [5Kb)) is a dark 
resonance that is adjacent to another dark resonance at 
Fg — ?> Fe — 2 (Fig. inia)), and hence its contrast is 
greater than for the resonance at the Fg = 2 ^ F^ = 2 
transition (Fig. [5Kc)), which is a dark resonance that is 
next to a resonance at Fg = 2 —^ F^ — 3, which is ex- 
pected to be bright. Figure SKd) shows a bright reso- 
nance, which is also expected because Fg < F^. The 
contrast of the resonance in Fig. [3Kd) is about 0.15%, 
but the agreement between experiment and calculation 
is still quite good. Figure igd) describes the *^^Rb Fg = 
2 — > Fe = 3 transition, which should also be a bright res- 
onance because Fg < Fe- However, the measured curve, 
which was recorded at the laser frequency that corre- 
sponded to maximum LIF, shows a dark resonance. We 
can understand why by studying Figure [51 

Figure [nKa)-(d) shows resonances obtained when the 
laser detuning with respect to the Fg = 2 ^ F^ = 3 
transition was -260 MHz, -160 MHz, -60 MHz, and 
-1-240 MHz, respectively. A detuning of -360 MHz with 
respect to the Fg = 2 ^ Fe = 3 transition would cor- 
respond exactly to the Fg — 2 ^ F^ — 2 transition, 
which is expected to be dark and is shown in Figure [SKc). 
The dark resonance at Fg = 2 — > Fe = 2 transition 
is clearly much stronger than the bright resonance at 
the Fg = 2 — > Fe = 3 transition, and overwhelms the 
bright resonance when the laser is tuned exactly to the 
Fg = 2 ^ Fe = 3 transition, or to the maximum in LIF 
near that transition at a slightly lower frequency. How- 
ever, as can be seen in Fig. [SI when the laser is tuned past 
the Fg = 2 — > Fe = 3 transition, further away from the 
Fg — 2 ^ Fe — 2 transition, a bright resonance can be 
observed. The respective theoretical calculations show 
the same behavior with good agreement. As far as we 
know such a detailed comparison of calculated and mea- 
sured non-linear magneto-optical resonance profiles as a 
function of laser detuning has not been presented previ- 
ously, and it is satisfying to note that a careful theoretical 
treatment can reproduce well the resonance profiles at an 
unresolved transition for laser excitation frequencies that 
differ by only 100 MHz. 

We also studied how the resonance shape varies as a 
function of laser power density. Figure [7| shows resonance 
signals for ®^Rb when the Fg = 3 — > Fe = 2 transition 
is excited at different laser power densities. The depen- 
dence of contrast on laser power density for all transitions 
can be extracted from plots similar to the ones in Fig. [3 
and presented as in Fig. [51 which shows the resonance 
contrast as a function of laser power density for the dark 
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FIG. 4: Fluorescence intensity versus magnetic field for Rb at Di excitation. Filled squares, experiment; solid line, theory. 
Note that the vertical axis is identical for (a)-(c), but markedly differs in (d). The excited state, total angular momentum of 
the ground Fg and excited states Fe of the transition, laser power density 7, and laser beam cross-section (FWHM) 5 are given 
in each panel. 



resonances of (a)*^Rb and (b)*''Rb and (c) for the bright 
resonance at the Fg = 1 ^ Fe — 2 transition of ^^Rb. 
The contrast C is defined as C = \Ires — Ifar\/Ifar, where 
Ires is the LIF intensity at the resonance at zero magnetic 
field and I far is the LIF intensity far from the resonance, 
where the curve of LIF versus magnetic field becomes 
fiat or has an extremum. Markers represent experimen- 
tally measured points, and fines represent the results of 
theoretical calculations. It can be seen that the model 
predicts the resonance contrast very well over more than 
two orders of magnitude in laser power density, and even 
at laser power densities at which one of the model's as- 
sumptions is not expected to be fully adequate, namely, 
the validity of describing the transit relaxation by a single 
rate constant [ssj . 

Figure [H shows the results of a study of how the reso- 
nance contrast and full width at half maximum (FWHM) 



depend on the transit relaxation rate. Transit relaxation 
is the process whereby atoms relax by leaving the laser 
beam. Thus, the transit relaxation rate is proportional 
to the inverse of the beam diameter. In Fig. [9l we plot 
the resonance contrast and full width at half maximum 
(FWHM) as a function of the inverse square root of the 
cross-sectional area of the beam, which is essentially the 
inverse of the beam diameter. The results are for the 
Fg = 2 — > Fg = 1 transition of ^^Rb. Again markers 
show the results of experimental measurements, whereas 
the line shows the result of calculations. The resonance 
width increases monotonically with relaxation rate, and 
the calculations predict the resonance width reasonably 
well over the range of relaxation rates studied. 
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FIG. 5: Fluorescence intensity versus magnetic field for Rb at Di excitation. Filled squares, experiment; solid line, theory. 
Note that the vertical axis is identical for (a)-(c), but markedly differs in (d). The excited state, total angular momentum of 
the ground Fg and excited states Fe of the transition, laser power density /, and laser beam cross-section (FWHM) S axe given 
in each panel. 



V. CONCLUSION 

The bright and dark magneto-optical resonances Di 
transition of natural rubidium have been carefully mea- 
sured with special efforts to eliminate systematic errors, 
and a detailed theoretical treatment has been applied 
to describe the experimental results in a system that is 
only partially resolved under Doppler broadening. The 
model wa.s based on the optical Bloch equations for the 
density matrix and included averaging over the Doppler 
profile, mixing of magnetic sublevels in an external mag- 
netic field, participation of all nearby transitions, and 
a detailed treatment of the coherence properties of the 
exciting laser radiation. The calculation based on our 
model accurately reproduces the experimental signals at 
each hyperfine transition and over a wide range of laser 
power densities and beam diameters. A bright resonance 
was observed at the Fg = 2 ^ = 3 transition of ^^Rb, 



but displaced from the exact position of the transition 
due to the influence of the nearby Fg = 2 Fe — 2 
transition, which is a dark resonance whose contrast is 
almost two orders of magnitude larger than the contrast 
of the bright resonance at the Fg ^ 2 ^ Fe = 3 transi- 
tion. The amplitude of the bright resonance is so small 
because the bright resonance is formed by a bright state 
that is strongly and continuously coupled to the laser field 
in a partially open transition, and as a result, this popu- 
lation of atoms in this state is heavily depleted. In con- 
trast, the dark resonance is formed by a dark state which 
is completely decoupled from the light field and which 
can efficiently maintain its population even for an open 
transition. Even in this very delicate situation, the the- 
oretical model described the experimental signals quite 
well at different laser detunings. As far as we know, the 
the bright resonance expected at the Fg = 1 ^ Fe = 2 
transition of *^Rb had not been observed before. 
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FIG. 6: Fluorescence intensity versus magnetic field for tlie Fg — 2 ^ Fe — 3 transition of ®^Rb at various values of the laser 
detuning with respect to the exact position of the Fg — 2 ^ Fe — 3 transition. Filled squares, experiment; solid line, theory. 



We believe that this study tests the theoretical model 
based on the optical Bloch equations to its fullest poten- 
tial because experiment and theory were compared for all 
eight hyperfine transitions of two isotopes at various laser 
power densities, laser beam diameters and laser frequen- 
cies. Agreement between experiment and theory over 
this wide range of situations and experimental conditions 
was excellent. Previous studies that presented compar- 
isons between experiment and theory for magneto-optical 
resonances in similar atomic systems have not had to 
take into account the Doppler effect because they used 
atomic beams [ill, [13, worked in systems resolved under 
Doppler broadening [22| , have not presented comparisons 
of theory and experiment over such a wide range of situ- 
ations [131 , or have not included in the theoretical model 
a detailed treatment of the coherence properties of the 



laser radiation or effects such as mi xing of magnetic sub- 
levels in an external magnetic field [l7| . 
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